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MAGNETIC RESONANCE 
IMAGING WITH NESTED 
GRADIENT PULSES 

Background of Invention 

[0001] The field of the invention is nuclear magnetic resonance imaging methods and 
systems. More particularly, the invention relates to a method and system for nesting 
gradient pulses to allow increased gradient slew rates and reduced peripheral nerve 
stimulation. It will be appreciated, however, that the-invention is also amenable to 
other like applications. 

[0002] Magnetic resonance imaging is a diagnostic imaging modality that does not rely 
on ionizing radiation. Instead, it uses strong (ideally) static magnetic fields, radio- 
frequency (RF) pulses of energy and magnetic field gradient waveforms. More 
specifically, MR imaging is a non-invasive procedure that uses nuclear magnetization 
and radio waves for producing internal pictures of a subject. Three-dimensional 
diagnostic image data is acquired for respective "slices" of an area of the subject 
under investigation. These slices of data typically provide structural detail having a 
resolution of one (1) millimeter or better. 



[0003] 



When utilizing NMR to produce images, a technique is employed to obtain NMR 

signals from specific locations in the subject. Typically, the region that is to be imaged 

(region of interest) is scanned by a sequence of NMR measurement cycles, which vary 

according to the particular localization method being used. The resulting set of 

received NMR signals are digitized and processed to reconstruct the image using one 

of many well-known reconstruction techniques. To perform such a scan, it is, of 

course, necessary to elicit NMR signals from specific locations in the subject. This is 

accomplished by employing gradient magnetic fields denoted (G , G 

SLICE PHASE 
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ENCODE ' and G READOUT ^ which have the same direction as the polarizing field B , 
but which are configured as needed to select the slice, phase encode and readout to 
facilitate the imaging. By controlling the strength of these gradients during each NMR 
cycle, the spatial distribution of spin excitation can be controlled and the location of 
the resulting NMR signals can be identified. 

[0004] The magnetic field gradient subsystem of an MR! system is perhaps the most 
critical subsystem in defining the utility of a scanner. A more powerful gradient 
subsystem, in general, yields greater applications capability. Gradient field 
performance is dependent on two factors: maximum gradient amplitude and gradient 
pulse slew-rate, Gradient amplitude is determined by the current that the gradient 
amplifiers produce in the gradient coils, and gradient slew rate is the rate at which the 
gradient amplifiers can change the gradient amplitude. 

[0005] In many circumstances, the only factor of importance in the generation of a 
gradient field pulse is the integral of gradient amplitude over the duration of the 
gradient pulse (i.e. the gradient pulse area). This is true, for example, with slice-select 
refocusing, phase-encoding, velocity or flow compensation, spoiling, rewinding and 
readout defocusing gradient pulses. Since the shortest duration gradient pulse of a 
given area provides the greatest flexibility in selecting pulse sequence echo time (TE) 
and pulse sequence repetition time (TR), it is highly desirable for the MRI system to 
produce these gradient pulses with the minimum pulse duration possible given the 
prescribed pulse area. 

[0006] 

For a given gradient subsystem, the minimum pulse duration is obtained with a 
triangular pulse. In such a case, a gradient amplifier is changing the gradient 
amplitude rapidly until the amplitude reaches a predetermined value. The gradient 
amplifier then returns the amplitude to zero at a rapid rate. The maximum rate that 
the gradient amplifier can change the gradient amplitude is established by the slew 
rate of the gradient amplifier given in units of Gauss per centimeter per millisecond 
(G/cm/ms). Ideally, a triangular pulse is the shortest possible gradient pulse for a 
given gradient amplifier since the gradient amplifier is changing the gradient 
amplitude at the maximum rate, established/limited by design. A real gradient 
amplifier, however, has a limit to the amplitude of gradient pulses it can produce. 
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Consequently, when the gradient amplifier must formulate a gradient pulse having an 

area greater than that of a triangular pulse with the maximum allowable amplitude, 

the gradient pulse providing optimal duration becomes a trapezoid. As with the 

triangular pulse the gradient amplitude rises at a rapid rate until the maximum 

gradient, G , is reached. The gradient amplifier then provides constant gradient 
max 

amplitude for a period PW, followed by a rapid return to zero amplitude. For the 
minimum duration gradient pulse, the constant amplitude is applied at the design 
limit of the gradient amplifier. 

[0007] The desired area of a gradient pulse is typically determined by scan parameters 
such as the field-of-view and slice thickness, which are input by the operator just 
prior to a scan. Once the desired area is known, then the MRI system uses the slew 
rate and maximum amplitude of the gradient amplifier to determine the timing and 
amplitude of the gradient pulse. If the is lew rate and maximum amplitude of the 
gradient system are sufficiently high, a triangular gradient waveform, will be 
produced. Otherwise, a longer trapezoidal gradient waveform must be produced with 
a consequent lengthening of the minimum possible pulse sequence echo time TE and 
pulse sequence repetition time (TR). 

[0008] While one could increase the power of the gradient amplifiers such that the 

maximum prescribed gradient area can be produced with a triangular gradient pulse 
regardless of its size, this may not solve the problem, and may present others. 
Namely, due to physiological effects on the patient, because constraints are placed on 
maximum gradient switching speed (slew rate) for the gradient fields allowed. Time- 
varying magnetic fields induce currents in conductive materials and rapidly changing 
magnetic field gradients can induce currents in a patient being imaged. Under some 
circumstances, these induced currents can stimulate nerves. Therefore, every MRI 
pulse sequence employed for human patients must conform to one or more magnetic 
field rate of change limitations in accordance with FDA regulations. Thus, it is now 
evident, that even though the gradient amplifiers can deliver a desired triangular 
gradient pulse of prescribed area to meet the desired requirements, physiological 
limits may preclude its use. Current MRI systems, therefore, assume the worst 
possible circumstances and limit the gradient slew rate accordingly. 
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[0009] Most physiological limits placed on the gradient field rate of change are not a 

single fixed value. Instead, the limit changes as a function of the "transition time" (i.e. 
the time interval over which the change in gradient field occurs). The reason for 
allowing higher rates of change (i.e. dB/dt) as the transition time decreases is related 
to the fact that the electrical sensitivity of neurons decrease with increasing frequency. 
J. P. Reilly of the Johns Hopkins University Applied Physics Lab has modeled the 
response of nerve cells and produced an equation predicting the dB/dt threshold for 
peripheral (PNST) and cardiac nerve stimulation as a function of dB/dt and pulse 
duration. It should be noted that cardiac stimulation occurs at dB/dt levels about 10 
times that of PNST, therefore, a wide margin of safety is realized. The Reilly PNST 
equation, known as the "Reilly Curve," is the basis for the FDA physiologic limits on 
dB/dt. 

[0010] All gradient coil designs intended for human use will have a physiologic limit 

given by the Reilly equation. The slew rate, which gives the limit, however, will depend 
on the effective length of the coil. The physiologic slew rate limit is determined by 
dividing the Reilly limit by the effective length of the gradient coil, L Note that, the 
effective length L is not necessarily the true length of the coil. Rather, the effective 
length is the distance between the iso-center of the gradient coil and the point within 
the coil that exhibits the maximum dB/dt exposure for the patient (or operator). This 
distance should be corrected for deviations from non-linearity. For a cylindrical 
gradient coil, this distance is roughly equal to the distance from the iso-center of the 
gradient coil to the location of maximum field variation caused by the coil in the Z 
direction, in the radial dimensions, however, this distance is equal to the patient-bore 
radius because the location of maximum field variation is within the walls of the MR 
system. 

[001 1] Most NMR scans currently used to produce medical images require many minutes 
to acquire the necessary data. The reduction of this scan time is an important 
consideration, since reduced scan time increases patient throughput, improves patient 
comfort, and improves image quality. Thus, it can be seen that there is a need to 
provide minimum duration gradient pulses that do not operate beyond the 
physiological limits. 
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Summary of Invention 

[0012] The above discussed and other drawbacks and deficiencies are overcome or 
alleviated by a system and method for decreasing gradient field pulse sequence 
duration and reducing peripheral nerve stimulation with known gradient pulse areas 
for a magnetic resonance imaging system. The method comprising: receiving a first 
desired area corresponding to a first pulse; obtaining a second desired area 
corresponding to a second pulse; selecting the first pulse as a nested pulse if the first 
desired area is smaller than the second desired area, and establishing the second 
pulse as a nesting pulse, otherwise selecting the second pulse as the nested pulse and 
establishing the first pulse as the nesting pulse. The method also includes: 
determining an amplitude and pulse duration for the nested pulse and ascertaining an 
amplitude and pulse duration for the nesting pulse. Finally, the method includes 
arranging a plurality of gradient field pulse sequences to include the nested pulse and 
the nesting pulse. 

[001 3] Also disclosed is a storage medium encoded with a machine-readable computer 
program code including: instructions for causing a computer to implement the 
abovementioned method for decreasing gradient field pulse sequence duration and 
reducing peripheral nerve stimulation with known gradient pulse areas for a magnetic 
resonance imaging system. 

[001 4] Further disclosed is a computer data signal comprising code configured to cause a 
processor to implement the abovementioned method for decreasing gradient field 
pulse sequence duration and reducing peripheral nerve stimulation with known 
gradient pulse areas for a magnetic resonance imaging system. 

[001 5] The above discussed and other features and advantages of the present invention 
will be appreciated and understood by those skilled in the art from the following 
detailed description and drawings. 

Brief Description of Drawings 

[001 6] Referring to the exemplary drawings wherein like elements are numbered alike in 
the several Figures: 

[001 7] FIG. 1 depicts an exemplary MRI system; 
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[0018] FIG. 2 is a graphic representation of a conventional gradient sequence; 

[0019] FIG. 3 is a graphic representation of a gradient sequence employing an exemplary 
embodiment; 

[0020] FIG. 4 is a graphic representation of a gradient sequence employing an alternative 
embodiment; 

[0021] FIG. 5 is a graphic representation of a gradient sequence employing yet another 
alternative embodiment; 

[0022] FIG. 6 is a flowchart depicting the method of acquisition for an exemplary 
embodiment; 

[0023] FIG. 7 depicts a triangular pulse area; 

[0024] FIG. 8 depicts a trapezoidal pulse area; 

[0025] FIG. 9 depicts a nested triangular pulse and a trapezoidal nesting pulse; and 

[0026] FIG. 1 0 is a graphic representation of a gradient sequence employing an 
£H exemplary embodiment. 

Detailed Description 

[0027] Disclosed herein is another method and system for reducing gradient pulse 

durations and thereby reducing TE and TR. The method includes nesting of selecting 
gradient pulses to reduce the overall length of the set of gradient pulses and thus the 
reducing the total cycle time for acquiring the images. 



[0028] 



To improve MR methods and scans, it is desirable to minimize the duration of the 
gradient pulses. Increasing the slew rate will have the effect of reducing pulse 
duration. However, important physiological limits must be considered. Slew rates that 
create dB/dt above threshold levels, established for example by the Reilly equation, 
may induce nerve stimulation in patients. For any gradient coil, dB/dt is directly 
proportional to the slew rate, and varies as a function of position in and around the 
gradient coil. The geometry of the coil will also have a substantial effect. For example, 
the maximum dB/dt in a shorter cylindrical coil will be lower than the maximum dB/dt 
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in a longer coil. Moreover, there is an upper limit to the amount of usable electrical 
current (and by extension power) for all gradient subsystems, and this upper limit is 
defined solely by the threshold of physiological stimulation and the maximum 
gradient pulse area for all applications. 

[0029] A method for calculating the shortest possible gradient pulse for a prescribed 

gradient pulse is described in U.S. Patent number 6,1 98,282. In this patent, a method 
is described to calculate the slew rate and am plitude of each gradient pulse, and to 
minimize the gradient pulse duration while addressing physiologic effects. The 
methods described therein, enables the gradient amplitude to be determined to 
achieve the minimum duration gradient pulse with the given slew-rate limitation. 
Thereby the echo time (TE) and repetition time (TR) can be optimized for the specified 
pulse sequence. Moreover, the gradients produced by the gradient subsystem were 
de-rated to avoid operation beyond the Reilly curve, yet provide minimum duration 
gradient pulses. The approach disclosed herein, improves upon the abovementioned 
reference by nesting selected portions of the gradient waveforms, and thereby further 
reduces gradient pulse sequence durations and thus, further reducing TE and TR. The 
method includes nesting of selected gradient slews and pulses to reduce the overall 
duration of the gradient pulse sequences and thus the reducing the total cycle time 
for acquiring the images. 

[0030] Referring to FIG. 1 , there is shown the major components of an exemplary MRI 

system, within which an exemplary embodiment may be implemented. The operation 
of the system is controlled from an operator console 100, which includes a keyboard 
and control panel 1 02 and a display 1 04. The console 1 00 communicates through a 
link 1 1 6 with a separate computer system 1 07 that enables an operator to control the 
production and display of images on the screen 1 04. The computer system 1 07 
includes a number of modules, which communicate with each other through a 
backplane. These include an image processor module 1 06, a CPU module 1 08 and a 
memory module 113, known in the art as a frame buffer for storing image data arrays. 
The computer system 1 07 is linked to storage media 111 and 1 1 2, depicted as disk 
storage and a tape drive respectively for storage of image data and programs, and it 
communicates with a separate system control 122 through a high speed serial link 
115. 
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[003 1 ] The system control 1 22 includes a set of modules connected together by a 

backplane 1 1 8. These include a CPU module 1 1 9 and a pulse generator module 121, 
which connects to the operator console 1 00 through a serial link 125. It is through 
this link 125 that the system control 1 22 receives commands from the operator that 
indicate the scan sequence that is to be performed. As will be described in more detail 
below, the operator enters parameters, which indicate the prescribed scan. From these 
parameters, a pulse sequence is calculated and downloaded to the pulse generator 
module 121. 

[0032] The pulse generator module 121 operates the system components to carry out the 
desired scan sequence. It produces data, which indicates the timing, strength and 
shape of the RF pulses that are to be produced, and the timing of and length of the 
data acquisition window. The pulse generator module 121 connects to a set of 
gradient amplifiers 127, to indicate the timing and shape of the gradient pulses to be 
produced during the scan. The pulse generator module 121 also receives patient data 
from a physiological acquisition controller 1 29 that receives signals from a number of 
different sensors connected to the patient, such as ECG signals from electrodes or 
respiratory signals from a bellows. Finally, the pulse generator module 121 connects 
to a scan room interface circuit 1 33, which receives signals from various sensors 
associated with the condition of the patient and the magnet system. It is also through 
the scan room interface circuit 1 33 that a patient positioning system 1 34 receives 
commands to move the patient to the desired position for the scan. 

[0033] The grac |j ent waveforms produced by the pulse generator module 121 are applied 

to a gradient amplifier 1 27 comprised of G , G , and G amplifiers. Each gradient 

x y z 

amplifier 127 excites a corresponding gradient coil in an assembly generally 
designated 139 to produce the magnetic field gradients used for position encoding 
acquired signals. The gradient coil assembly 1 39 forms part of a magnet assembly 
141, which includes a polarizing magnet 140 and a whole-body RF coil 152. A 
transceiver module 1 50 in the system control 1 22 produces pulses, which are 
amplified by an RF amplifier 151 and coupled to the RF coil 1 52 by a transmit/receive 
switch 1 54. The resulting signals radiated by the excited nuclei in the patient may be 
sensed by the same RF coil 1 52 and coupled through the transmit/receive switch 1 54 
to a preamplifier 1 53. The amplified MR signals are demodulated, filtered, and 
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digitized in the receiver section of the transceiver 1 50. The transmit/receive switch 
1 54 is controlled by a signal from the pulse generator module 1 21 to electrically 
connect the RF amplifier 1 51 to the RF coil 1 52 during a transmit mode and to 
connect the preamplifier 1 53 during a receive mode. The transmit/receive switch 1 54 
also enables a separate RF coil 1 52 (for example, a head coil or surface coil) to be 
used in either the transmit mode or receive mode. 

[0034] The MR signals picked up by the RF coil 1 52 are digitized by the transceiver 
module 1 50 and transferred to a memory module 1 60 in the system control 1 22. 
When the scan is completed and an entire array of data has been acquired in the 
memory module 1 60, an array processor 161 operates to Fourier transform the data 
into an array of image data. This image data is conveyed through the serial link 1 1 5 to 
the computer system 1 07 where it is stored in a storage medium 111 or 1 1 2 such as 
disk memory or tape drive. The storage medium 1 1 1 and 1 1 2 could be various storage 
methodologies, such as disk, static memory, solid state, removable media, and the 
like, as well as combinations including at least one of the foregoing. In response to 
commands received from the operator console 100, this image data may be archived 
on the tape drive, or it may be further processed by the image processor 1 06, and 
conveyed to the operator console 1 00 and presented on the display 1 04. 

[0035] Referring still to FIG. 1 the NMR signal produced by the subject is picked up by the 
receiver coil 1 52 and applied through the preamplifier 1 53 to the input of a 
transceiver 1 50. The received signal is at or around the Larmor frequency, and this 
high frequency signal is down converted in a two-step process, which first mixes the 
NMR signal with a carrier signal and then mixes the resulting difference signal with a 
reference signal. The down converted NMR signal is applied to the input of an analog- 
to-digital (A/D) converter, which samples and digitizes the analog signal and applies 
it to a digital detector and signal processor which produces in-phase (I) values and 
quadrature (Q) values corresponding to the received NMR signal/The resulting stream 
of digitized I and Q values of the received signal are output through backplane 1 1 8 to 
the memory module 1 60 and array processor 1 61 where they are employed to 
reconstruct an image. 

[0036] 

Referring to FIG. 2, a graphic representation of an exemplary gradient sequences 
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are depicted. The gradient magnetic field sequences denoted (G , G , and 

SLICE PHASE 

G READ0U j ) which have the same direction as the polarizing field B , but which 

have a gradient along the respective z, y and x axes. It will be appreciated, that G 

.. , G , and G may be applied on any physical gradient axes G , C 

slice readout phase x 

, G „ or any linear combination thereof, depending on the selected orientation for 
y z 

the slice. The slice select pulse sequence, G , 20 includes, but is not limited to, a 

slice 

slice select pulse 22, a rephaser pulse 24, and a killer or dephaser pulse 26. The 

phase encoding pulse sequence, G , 30 includes, but is not limited to, the phase 

phase 

encoding pulse(s) 32 and phase rewinder pulse(s) 34. Lastly, the readout pulse 

sequence, G 40 includes, but is not limited to, a prephaser pulse 42, a 

readout 

readout pulse 44, and a killer or dephaser pulse 46. 



,0 [0037] 

ff' Referring now to FIG. 3, a graphical representation of a gradient sequence 

UJ 

employing an exemplary embodiment of the invention is depicted. The gradient 
| magnetic field pulses as described and depicted from FIG. 2 are reconfigured in a 

manner, which ensures that the pulses for slice selection gradient field and the 
readout field are nested or coincident. More particularly, the slewing portions of the 
pulses of the gradient field sequences are arranged to minimize any occurrence of 
p more than one gradient field slewing simultaneously. In an exemplary embodiment, 

the pulse sequences e.g., 20, 30, and 40 are initiated with the complete slice select 
pulse 22 in the slice select pulse sequence 20. Thereafter, the initial (e.g., first in time 
and denoted generally as 52 in FIG. 9), slew 52 of the slice select rephaser pulse 24 is 
generated, with the readout prephaser pulse 42, nested in the non-slewing e.g., flat 
top or unchanging field, portion of the slice select rephaser pulse 24. Thereafter, the 
initial slew 52 for the readout pulse 44 of the readout pulse sequence 40 starts just 
after the completion of the terminal slew (e.g., later in time and denoted generally as 
54 in FIG 9) of the slice select rephaser pulse 24. This embodiment places a the entire 
phase encoding pulse(s) 32 and phase encoding rewinder pulse(s) 34 of the phase 
pulse sequence 30 within the readout pulse 44, during which the echo is being 
captured. The phase encoding pulse(s) 32 and their corresponding rewinder pulse(s) 
34 all occur during the non-slewing portion of the readout pulse 44. At the 
termination of the readout pulse 44 (following echo acquisition), the slice select 
dephaser pulse 26 of the slice select pulse sequence 20 slews with the readout 
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dephaser pulse 46 of the readout pulse sequence 40 nested therein. This embodiment 
is beneficial because it represents the simplest configuration for the gradient pulse 
sequences e.g., 20, 30, and 40. Furthermore, it will be appreciated that the pulse 
sequences could readily be rearranged with the slice select rephaser pulse 24 and 
dephaser pulse 26 nested within the non-slewing portion of the readout prephaser 
pulse 32 and readout dephaser pulse 46 respectively. 

[0038] FIG. 4 depicts a graphical representation of a gradient sequences e.g., 20, 30, and 
40 for an alternative embodiment. It should be noted that this sequence is similar to 
that described above, except that the slice select rephaser pulse 24 and dephaser 
pulse 26 are nested during the non-slewing portions of the readout prephaser pulse 
42 and dephaser pulse 46 respectively. Additionally, the initial slew 52 for the phase 
encoding pulse(s) 32 and the terminal slew 54 for the phase rewinder pulse(s) 34 are 
placed outside the readout pulse 44 of the readout pulse sequence 40. Namely, the 
initial slew 52 for the phase encoding pulse(s) 32 is placed in between the terminal 
slew 54 of the readout prephaser pulse 42 and the initial slew 52 of the readout pulse 
44. In this embodiment, the gradient magnetic field pulses as described and depicted 
from FIG. 2 are once again, reconfigured in a manner, which ensures that the pulse for 
slice selection gradient field and the readout field are nested or coincident. Once 
again, in this alternative embodiment, the sequences are initiated with the complete 
slice select pulse 22. Thereafter, the initial slew of the readout prephaser pulse 42 is 
generated with the slice select rephaser pulse 24 nested in the non-slewing portion of 
the readout prephaser pulse 42. Thereafter, the initial slew 52 for the phase encoding 
pulse(s) 32 starts just after the completion of the readout prephaser pulse 42. During 
the non-slewing portions of the phase encoding pulse(s) 32 and their corresponding 
rewinder pulse(s) 34, the initial slew 52 and terminal slew 54 portions, respectively of 
the readout pulse 44 are generated. At the termination of the phase rewinder pulse(s) 
34, the readout dephaser pulse 46 slews with the slice select dephaser pulse 26 
nested therein. 

[0039] Once again, it will be appreciated that the sequences described above could 
readily be rearranged with the readout prephaser pulse 42 and dephaser pulse 46 
nested within the non-slewing portion of the slice select prephaser pulse 24 and slice 
select dephaser pulse 26 respectively. 
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[0040] FIG. 5 depicts a graphical representation of a gradient sequence for yet another 
alternative embodiment. Once again, this sequence is similar to that described above, 
except that the slewing portions of the phase encoding pulse(s) 32 and phase 
rewinder pulse(s) are displaced even more from the readout pulse 44. In particular, 
The initial slew 52 of the phase encoding pulse(s) 32 is nested within both the slice 
select rephaser pulse 24 and readout prephaser pulse 42, while the terminal slew 54 
of the phase encoding rewinder pulse(s) 34 is nested during the non-slewing portions 
of the slice select dephaser pulse 26 and readout dephaser pulse 46 respectively. 
Once again, in this alternative embodiment, the sequences are initiated with the 
complete slice select pulse 22. Thereafter, the initial slew 52 of the slice select 
rephaser pulse 24 is generated with the complete readout prephaser pulse 42 nested 
in the non-slewing portion of the slice select rephaser pulse 24. The slice select 
rephaser pulse 24 non-slewing portion is extended for the duration needed to 
complete the readout prephaser pulse 42. Unique in this embodiment, the initial slew 
52 for the phase encoding pulse(s) 32 starts just after the completion of the initial 
slew 52 of readout prephaser pulse 42. The readout prephaser pulse 42 is limited to 
not slewing during the duration of the initial slew 52 of the phase encoding pulse 32. 
During the non-slewing portions of the phase encoding pulse(s) 32 the terminal 
slewing portions 54 of the slice select rephaser pulse 24 and the initial slew 52 of the 
readout pulse 44 are generated. Immediately following the completion of the initial 
slew 52 of the readout pulse 44, the terminal slew 54 phase encoding pulse(s) 34 is 
completed. Likewise, at the end of the echo acquisition/the initial slew 52 of the 
phase encoding rewinder pulse 34 is generated, followed thereafter by the terminal 
slew 54 of the readout pulse 44, and then the initial slew 52 of the slice select 
dephaser pulse 26. Nested completely within the non slewing portion of the slice 
select dephaser pulse 26 is the complete readout dephaser pulse 46. Moreover, 
nested in the non-slewing portion of the readout dephaser pulse 46 is the terminal 
slew 54of the phase encoding rewinder pulse(s) 34. Lastly, at the completion of the 
terminal slew 54 of the readout dephaser pulse 46, the terminal slew 54 of the slice 
select dephaser pulse 26 is generated. 

[0041] 

It should be appreciated that yet another embodiment could be conceived as a 
variation on the embodiment described above and depicted in FIG. 5. In this 
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alternative embodiment, the slice select rephaser pulse 24 may instead be nested 
within the duration of the readout prephaser pulse 42, with the initial slew 52 of the 
phase encoding pulse(s) 32 being generated during the non-slewing portions of the 
slice select rephaser pulse 24. Likewise, the slice select dephaser pulse 26 may 
instead be nested within the duration of the readout dephaser pulse 46 with the 
terminal slew 54 of the phase encoding rewinder pulse(s) 34 being generated during 
the non-slewing portions of the slice select dephaser pulse 26. 

[0042] Referring now to FIG. 6, a flow chart of the method employed in an exemplary 

embodiment is depicted. Disclosed herein is a method to configure a set of gradient 
field sequences including complete overlap of the phase encoding pulses with the 
readout window as depicted in FIG. 3. The method addresses a configuration for the 
slew rates, magnitudes and durations of the sequence pulses based upon particular 
constraints of the MRI system 10, the gradient amplifiers 1 27, as well as the 
limitations required by the Reilly equation to avoid peripheral nerve stimulation. For 
example, as stated earlier, the gradient amplifier 127 includes limitations as to the 
maximum field amplitudes and slew rates that it can produce. Additionally, to achieve 
desirable echo acquisition and for the reasons discussed above, it is known that for a 
given sequence the area of the various pulses is a primary consideration. 

[0043] it S | 10U | C | k e appreciated that the method disclosed herein does not address the 
formulation of the slice select pulse 22, the phase encoding pulse(s) 32, the phase 
rewinder pulse(s) 34 or the readout pulse 44. These pulses are selected as a function 
of numerous elements including, but not limited to, the type of gradient sequence to 
be formulated, the imaging to be performed, the characteristics of the MRI system 10, 
the gradient amplifiers 127 and the like. Therefore, it will be appreciated that the 
method disclosed herein is directed to the sequencing and shaping of the various 
gradient pulses based upon a known or computed desired area for each pulse. 
Further, it is well appreciated that the areas of the slice select rephaser pulse 24 and 
dephaser pulse 26 are about equivalent. Moreover, their combined area should 
preferably be about equivalent to the area of the slice select pulse 22. Similarly, the 
area of the phase encoding pulse(s) 32 should preferably be about equivalent to the 
area of the phase encoding rewinder pulse(s) 34. Further, it is well appreciated that 
the areas of the readout prephaser pulse 42 and dephaser pulse 46 are about 
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equivalent, and again, their combined area should preferably be about equivalent to 
the area of the readout pulse 44. 



[0044] Referring again to FIG. 6, a method to configure a set of gradient field sequences 
is described for a configuration for the slew rates, magnitudes and durations of the 
sequence pulses based upon particular constraints of the MRI system 1 0, the gradient 
amplifiers 127, as well as the limitations required by the Reilly equation to avoid 
peripheral nerve stimulation. For example: for a given MRI system 1 0 with the 
following constraints and known parameters: 

Known desired area for each pulse of the gradient field sequences 

Maximum field amplitude =A = 4 Gauss/centimeter (G/cm) 

Maximum slew rate - SR = nominally 1 50 milliTcsla/meter/second (mT/m/s) 

= (4 G/cm)/(268microseconds(jis)J 

= 0.0149 G/cm/us 



[0045] Minimum duration gradient pulse width for a known gradient amplifier is 4 y s in 
duration. Hence, a symmetric trapezoid pulse (or a minimum duration triangular pulse 
will actually have a flat top time of 8 [j s. Therefore, strictly speaking, all gradient 
pulses generated are trapezoidal because achieving a pure triangular pulse is beyond 
the capabilities of the MRI system 1 0. The derivation is presented addressing the 
triangular pulse because ideally, a triangular pulse is desired as it provides the 
maximum area within a minimum duration. 

[0046] 



Referring now to FIGs, 7 and 8, the area of a triangle is known to be: 
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areatnangje = 1 /£base*height = Vi 1^* Atnangic- 
where t ram p is the time duration of the ramp 
Atriangk is the amplitude of the triangle 

areamanguiar = 2C/2base*height) ^ 2(H t ramp *A lr iangie) = 2(16 t ramp *A lr i at igic) J 
area trianeilkr = t ran , p *A tr}0 „gi e - t ramp *A - (4 G/cm)/SR (G/(cm*us) * (4 G/cm) 
While the area of a rectangle is known to be: 

area re ctangle = tfiattop * A 

The desired area for a particular pulse is: areanns (G/cm * us) 
For a desired triangular pulse: 

If area D fcs <= (4 G/cm) / (0.0149 G/cm/us) * 4 G/cm, then use a triangular 
waveform. 

tramp = sqrt [area DES / SR] sqrt[(G us/cm)/(G/(cm*us))] 
AtrianeJc = sqrt LareaoEs (G us/cm) * SR (G/(cm*us)] 

For a desired trapezoidal pulse: 

If area > ((4 G/cm) / (0.0149 G/cm/us) + 8 us) * 4 G/cm, then use a trapezoidal 
waveform. 

Ain 3 ngic = A = 4 G/cm 

tramp = (4 G/cm) / SR (G/(cm*us) ) 

tfiattop = area D ES - (tramp + A triang i c )) 



[0047] Once again, it should be appreciated that in the embodiments presented herein, 
all pulses are treated as trapezoidal to address the constraints of the gradient 
amplifiers 127. Therefore, for those pulses that would otherwise be considered 
triangular, they are herein treated as trapezoidal with a flattop duration of 8 |j s. 

[0048] Referring once again to FIG. 6 and now to FIG. 9, in an exemplary embodiment 

exhibiting slice select rephaser/dephaser and readout rephaser/dephaser overlap the 
method includes a process to evaluate the desired area for each pulse. The desired 
areas for a first pulse e.g. slice select rephaser/dephaser 24, 26 and second pulse 
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e.g., readout rephaser/dephaser respectively are compared at process block 210 and 

decision block 220. The pulse exhibiting the smallest area is established as the nested 

pulse 60, while the pulse exhibiting the larger area is established as the nesting pulse 

70. Referring to process block 230, the amplitude and pulse durations for the nested 

pulse 60, whether triangular or trapezoidal, are determined. At process block 240, the 

amplitude and pulse durations of the nesting pulse 70 are determined. It will be 

appreciated that because of the gradient amplifiers 127 the nested pulse is 

trapezoidal despite even though, ideally, the nested pulse 60 should be configured as 

triangular. The flattop duration of the nested pulse 60 is therefore the minimum 

realizable pulse width of the gradient amplifier 127, in this instance, 8 microseconds. 

Moreover, it should also be appreciated that because of the fixed duration 

requirement for the nested pulse 60 during the non-slewing portion of the nesting 

pulse 70, the nesting pulse 70 is inherently trapezoidal. The flat top duration of the 

nesting pulse 70 is then set equal to the total pulse width of the nested pulse 60. For 

example: For the nested pulse 60, t , t . , and A , are known. For the 

ramp flattop triangle 

nesting pulse, 70 set the flattop or non-slewing duration to the total duration for the 
nested pulse 60; 

tfkuop_ncsting ~ Wied Pulse 

lflaLtop_ncsting = tramp_nested + tflallop nested + trump_neswd 

Note, t = 8 m s for a minimum duration nested pulse. Ideally, 0 u s for a 

flattop nested y 

triangular nested pulse 60. 

[0049] 



Process block 240 describes the process to compute the amplitude of the nesting 

pulse. To facilitate performing the computations, known quantities are the pulse area, 

t , and SR (slew rate), 

flattop 
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area D ES = arearectangie + arearnanguiar 
= A *t flattop + A*t ramp ) 

= A * (tramp + tflattop) 

=A*(A/SR + t llattop ) 

therefore, 

area DES * SR = A 2 + A*SR *lfi aU o P 
rearranging the equation, yields: 

A 2 - A*SR*t flattop - areaocs *SR = 0 
This result is a quadratic equation, which can readily be solved for A yielding: 

A = - SR*tft allot /2 +/- sqrt[SR*t na «op 2 -4 arca DH s *SR]/2 

Neglecting the incongruent solution, yields the amplitude, A 
Finally, to compute the ramp times: 

tramp ~ A / SR 

[0050] Returning now to FIG. 6, the process is completed by arranging the respective 
gradient pulse sequences e.g., 20, 30, and 40 as depicted at process block 250. 

[0051] Referring not to FIG. 1 0 is an exemplary depiction of the pulse sequences from 
FIG. 3 with illustrative pulse amplitudes and durations. This pulse sequence utilizes 
nested slice select rephaser and dephaser pulses 24, 26 and phase encoding pulse 32 
and rewinder pulse 34 overlap with the readout pulse to minimize dB/dt peripheral 
nerve stimulation effects. Slice select rephaser pulses 24 and dephaser pulses 26 are 
always triangular with pulse amplitudes of - 1 .91 366 G/cm with pulse ramp widths of 
132 m s to a minimum of +/- 0.146767 G/cm with pulse ramp widths of 12 M s. Phase 
encoding flat top widths are 8 ms in duration. Phase encoding pulses 32, 34 are 
always triangular with pulse amplitudes ranging from a maximum. of +/- 1.594516 
G/cm with pulse ramp widths of 1 08 p s. Phase encoding flat top widths are 8 p s in 
duration. Readout prephaser pulses 42 and dephaser pulses 46 are configured as 
trapezoidal with pulse amplitudes of - 1 .26025 G/cm with pulse ramp widths of 84 p 
s and a non-slewing, flat top duration of 272 p s. It should be noted, that for this 
exemplary sequence, the maximum slew rate of the gradient amplifier 127 is always 
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employed. Other slew rates and configurations of the sequences are possible. 

[0052] The system and methodology described in the numerous embodiments 

hereinbefore provides a system and method for nesting gradient pulses and thereby, 
increasing gradient slew rates and decreasing peripheral nerve stimulation. In 
addition, the disclosed invention may be embodied in the form of computer- 
implemented processes and apparatuses for practicing those processes. The present 
invention can also be embodied in the form of computer program code containing 
instructions embodied in tangible media, such as floppy diskettes, CD-ROMs, hard 
drives, or any other computer-readable storage medium such as 1 1 1 and 112, 
wherein, when the computer program code is loaded into and executed by a 
computer, the computer becomes an apparatus for practicing the invention. The 
present invention can also be embodied in the form of computer program code, for 
example, whether stored in a storage medium 111, 112, loaded into and/or executed 
by a computer, or as data signal 1 1 0 transmitted whether a modulated carrier wave or 
not, over some transmission medium, such as over electrical wiring or cabling, 
through fiber optics, or via electromagnetic radiation, wherein, when the computer 
program code is loaded into and executed by a computer, the computer becomes an 
apparatus for practicing the invention. When implemented on a general-purpose 
microprocessor, the computer program code segments configure the microprocessor 
to create specific logic circuits. 

[0053] While the invention has been described with reference to a preferred embodiment, 
it will be understood by those skilled in the art that various changes may be made and 
equivalents may be substituted for elements thereof without departing from the scope 
of the invention. In addition, many modifications may be made to adapt a particular 
situation or material to the teachings of the invention without departing from the 
essential scope thereof. Therefore, it is intended that the invention not be limited to 
the particular embodiment disclosed as the best mode contemplated for carrying out 
this invention, but that the invention will include all embodiments falling within the 
scope of the appended claims. Moreover, the use of the terms first, second, etc. do 
not denote any order or importance, but rather the terms first, second, etc. are used 
to distinguish one element from another. 
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